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Abstract 

We discuss relative roles played by the magnetic inversion symmetry breaking and the ferroelec- 
tric (FE) atomic displacements in the multiferroic state of YMnC-3. For these purposes we derive a 
realistic low-energy model, using results of first-principles electronic structure calculations and ex- 
perimental parameters of the crystal structure below and above the FE transition. Then, we solve 
this model in the mean-field Hartree-Fock approximation. We argue that the multiferroic state in 
YMn03 has a magnetic origin, and the centrosymmetric Pbnm structure is formally sufficient for 
explaining main details of the noncentrosymmetric magnetic ground state. The relativistic spin- 
orbit interaction lifts the degeneracy, caused by the frustration of isotropic exchange interactions in 
the ab plane, and stabilizes a twofold periodic noncollinear magnetic state, which is similar to the 
E-state apart from the spin canting. The noncentrosymmetric atomic displacements in the P2\nm 
phase reduce the spin canting, but do not change the symmetry of the magnetic state. The effect of 
the P2\nm distortion on the FE polarization AP a , parallel to the orthorhombic a axis, is twofold: 
(i) it gives rise to ionic contributions, associated with the oxygen and yttrium sites; (ii) it affects 
the electronic polarization, mainly through the change of the spin canting. The relatively small 
value of AP a , observed in the experiment, is caused by a partial cancelation of the electronic and 
ionic contributions, as well as different contributions in the ionic part, which takes place for the 
experimental P1\nm structure. The twofold periodic magnetic state competes with the fourfold 
periodic one and, even in the displaced P2\nm phase, these two states continue to coexist in a 
narrow energy range. Finally, we theoretically optimize the crystal structure. For these purposes 
we employ the LSDA+J7 approach and assume the collinear E-type antiferromagnetic alignment. 
Then, use the obtained structural information again as the input for the construction and solu- 
tion of the low-energy model. We have found that the agreement with the experimental data in 
this case is less satisfactory and |AP a | is largely overestimated. Although the magnetic structure 
can be formally tuned by varying the Coulomb repulsion U as a parameter, apparently LSDA+J7 
fails to reproduce some fine details of the experimental structure, and the cancelation of different 
contributions in AP a does not occur. 
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FIG. 1. (Color online) Schematic view on the orbital ordering and magnetic interactions in the 
ab plane of YMn03. and j£* denote nearest-neighbor magnetic interactions, associated with 
the compressed (<) and stretched (>) bonds (shown by the filled and hatched lines, respectively) 
in the ferroelectric P2\nm structure. Jf and are the second- neighbor interactions along the 
orthorhombic axes a and b, respectively. J3 is the strongest third-neighbor interaction. 

I. INTRODUCTION 

Multiferroic materials, where the ferroelectric (FE) activity coexists with some long-range 
antiferromagnetic (AFM) order, have attracted a great deal of attention.- A possibility to 
control magnetic properties by the electric field (and vice versa) makes these materials 
promising for the creation of a new generation of electronic devices. Moreover, there is a 
fundamental interest to the phenomenon of multiferroicity itself and the microscopic origin 
of coupling between magnetic and electric degrees of freedom. 

The orthorhombic rare-earth manganites i?Mn03 are one of the key compounds for un- 
derstanding the basic concepts of such coupling. There are two mechanisms of FE activity 
that are actively discussed today. One of them is the FE displacements in the noncen- 
trosymmetric twofold periodic E-type AFM (ttii) state (Fig. [T]), 2 which is caused by the 
exchange striction effects and leads to the alternation of the long and short Mn-Mn bonds 
in the ab plane of i?Mn03. This mechanism is nonrelativistic, and is believed to take place 
in the small- R compounds, such as HoMn03 and YMnCV Although the ttii arrangement 
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itself is sufficient for the magnetic inversion symmetry breaking, the FE displacements are 
believed to be important for stabilizing the E-state.- Another mechanism is associated with 
the cycloidal spin-spiral alignment, which couples to the lattice via the relativistic spin-orbit 
(SO) interaction.- The characteristic compound, whose properties are typically regarded in 
the context of this type of theories, is TbMnC>3. It forms nearly fourfold periodic magnetic 
texture.- 

In the recent publication (Ref.lf]), one of us has proposed that there is no conceptual differ- 
ence between twofold and fourfold periodic manganites, and it is not quite right to consider 
separately the relativistic and nonrelativistic mechanisms of the FE activity. Both types of 
magnetic structures are formed due to the complex competition, involving three main in- 
gredients: frustrated isotropic exchange interactions of the nonrelativistic origin, single-ion 
anisotropy, and Dzyaloshinskii-Moriya interactions. Moreover, the magnetic structure, re- 
sulting from this competition, appears to be more complex, where the collinear fN4- E-type 
AFM alignment as well as the cycloidal spin-spiral alignment can be regarded only as the 
first (and rather crude) approximation to the true magnetic ground state. Particularly, it 
was argued that, even in the fourfold periodic TbMnC>3, the FE activity is caused by the de- 
formation of the homogeneous spin-spiral alignment (rather than the spin-spiral alignment 
itself), and in this sense, the situation is similar to the twofold periodic HoMnC>3.- 

Regarding the twofold periodic magnetic ordering, realized in HoMnC>3 (which is similar 
to YMnC^), it was proposed that, although the FE distortion apparently play an impor- 
tant role in stabilizing this magnetic state, its symmetry is dictated by the competition of 
magnetic interactions in the centrosymmetric Phnm phase.- This competition is responsible 
for the magnetic inversion symmetry breaking, whereas the FE displacements should follow 
this lowering of the magnetic symmetry. 

There are two important developments, which we would like to address in this work. 
First, the experimental crystal structure, realized below and above the FE transition (which 
coincides with the transition to the E-state), has been reported recently for YMn03. 7 Par- 
ticularly, it was shown that the FE transition at around 30 K is accompanied by the lowering 
of symmetry from Pbnm till P2\nm. The experimental parameters of the crystal structure 
and the atomic positions have been reported for two characteristic temperatures, 21 K and 
50 K, correspondingly below and above the transition. Thus, following the idea of the pre- 
vious work (Ref. |o), we are able to clarify several aspects, related to the formation of the 



twofold periodic magnetic texture and the FE polarization in YMnC>3. Particularly, what is 
the role of the P2\nm atomic displacements in the formation of the twofold periodic mag- 
netic texture and the FE properties of YMnC>3? Second, we have performed the structural 
optimization for the E-type AFM state of YMnC>3, by using both LSDA and LSDA+[7 
techniques (where, LSDA stands for the local-spin-density approximation). It allows us 
to address another important problem: all previous attempts to theoretically optimize the 
crystal structure for the E-type AFM state systematically overestimated the value of the FE 
polarization by an order of magnitude.-^ On the other hand, the theoretical calculations, 
which use the experimental P2 1 nm structure, seem to reproduce also the observed polar- 
ization value.- Therefore, it is very important to understand the origin of this problem on a 
microscopic level. 

The rest of the paper is organized as follows. Sec. [11] is devoted to the analysis of the low- 
energy model, which was constructed by using the experimental parameters of the crystal 
structure. The method is briefly described in Sec. Ill Al Then, in Sec. Ill Bl we discuss the 
effect of the FE displacements in the P2\nm phase on interatomic magnetic interactions. 
The summary of results without relativistic SO coupling is given in Sec. Ill C\ where we 
investigate the effect of the FE displacement on the total energy and the electric polarization. 
In Sec. IIIDI we consider the effect of the relativistic SO interaction on the magnetic ground 
state, and in Sec. Ill El we discuss the value of the FE polarization. Then, in Sec. IHH 
we perform similar analysis, using results of structural optimization for the collinear E- 
type AFM state, and present a critical comparison with the results, obtained by using the 
experimental crystal structure. Finally, a brief summary of the work is given in Sec. IIVI 

II. CONSTRUCTION AND SOLUTION OF THE LOW-ENERGY MODEL 
A. Method 

For the analysis of electronic and magnetic properties of YMnOs, we employ the same 
strategy as in previous publications, devoted to multiferroic manganites.-^ First, we con- 
struct an effective low-energy model (more specifically - a multiorbital Hubbard model) 
for the Mn 3d bands and derive parameters of this model from the first-principles electronic 
structure calculations. For these purposes, we use the experimental parameters of the Pbnm 
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and P2\nm structure.- Then, we use this model for the analysis of the magnetic ground state 
and the FE polarization, which depends on the details of the magnetic and crystal structure 
of YMn0 3 . 

The model itself has the following form: 

It is formulated in the basis of Wannier orbitals for the Mn 3ci bands, by starting from 
the electronic structure in the local-density approximation (LDA).— In these notations, 
each Greek symbol stand for the combination of spin (s= f or \) and orbital (m= xy, yz, 
3z 2 —r 2 , zx, or x 2 —y 2 ) variables. The parameters of the crystal field and the SO interaction 
are included to the site-diagonal part of tf?. Off-diagonal matrix elements of if? with respect 
to the site indices (i and j) stand for the transfer integrals. U a p r s are the screened Coulomb 
interactions. The construction of the model (jTJ, definitions and details of calculations of 



the model parameters were explained in the review article (Ref. llOl ). The behavior of these 
parameters for the series of orthorhombic manganites was discussed in previous publications 
(Refs. |6| and|9|). The parameters, obtained for the P2\nm and Pbnm phases of YMn03, are 
collected in the supplemental materials.— 

After the construction, the model ([1]) is solved in the mean-field Hartree-Fock (HF) 
approximation:— 

{£(k) + v}|C nk )=^k|C nk ), 

where t(k) is the Fourier image of Uj and V is the HF potential. The phases of t(k) for the 
each pair of atoms in the primitive cell are specified in the following way: 

t TiTj (k) = 'x r •H /' k :T T (2) 

where (rj+Rj) is the position of the lattice site j in terms of its coordinate Tj within the 
primitive cell and the lattice translation Rj, and t(k) is the supermatrix: i(k) = \\i TiTj (k)\\. 
The electric polarization is calculated by using the Berry-phase formalism:— 

AP a = i E ZtAn, - ^ NiN ^ {Gal [7 a (oo) - 7b (0)] . (3) 

Here, the first term is the ionic contribution, which is caused by the displacements (At^q) of 
the atomic charges (Zi) away from the centrosymmetric positions. It is calculated within the 



unit cell with the volume V. The second term is the electronic contribution (in the following 
denoted as AP el ), which can be related to the difference of the Berry phases (7) in the 
process of the adiabatic symmetry lowering.— In the case of improper multiferroics, where 
the inversion symmetry is broken by some complex magnetic order, the 'adiabatic symmetry 
lowering' can be understood as the rotation of spins away from the centrosymmetric ferro- 
magnetic (FM) configuration. For the ionic term, the reference point was constructed by 
symmetrizing the atomic positions of the P2\nm structure, so that they would transform 
to each other according to the symmetry operations of the group Pbnm. In the process of 
such symmetrization, the unit-cell volume was kept constant. 

The Berry phases can be evaluated on a discrete grid of the k-points, generated by the 
N 1 xN 2 y<N 3 divisions of the reciprocal lattice vectors {G a } (0 < s a < N a ) 

k Sl so o„ = — Gi -|- —Go ~\~ ~ G^. 

1,2,3 Ni N 2 N 3 d 

This yields the construction of the type:— 

n 2 -i N3-1 

^ = -J2J2 ImlnP ^3, (4) 

s 2 =0 s 3 =0 

where similar expressions for 72 and 73 can be obtained by the cyclic permutation of the 
indices 1, 2, and 3. Then, f S2 s 3 can be related to the overlap matrices S = || (C n k|Cn'k')|| 
between occupied eigenvectors \C n k) in two neighboring k-points (for instance, k=k SliS2jS3 
and k'=k sl+ i !S2)S3 for 7^ etc.). If |C„k) is periodic in the reciprocal space, 

|Cnk+G a ) = |Cnk) (5) 



(the so-called periodic gauge), P S2 s 3 is given by the following expression™ 

2Vi-l 

P s 2 s 3 — J J d e t5 , (k sljS2jS3 , k sl _|_i jS2jS3 ), (6) 

si=0 

where, according to Eq. (JSJ), |C n k) i n the point k^ 1)S2ja3 can be replaced by the one in the 
point ko iS2i s 3 . However, if the kinetic part of the model Hamiltonian is given by Eq. (j5J), the 
eigenvectors are not periodic. Therefore, we should modify f S2 s 3 by including the additional 
multiplier: 

F S2S3 ->■ P S2S3 detS , (k A r 1)S2jS3 , k 0)S2jS3 ). (7) 

It describes the change of the phase at the zone boundary, which appears when we re- 
duce the k-space integration to the first Brillouin zone. Furthermore, since Eq. (jSJ) is no 



longer applicable, the eigenvectors \C n k) in the point k.Ni,s a ,ss should be obtained from in- 
dependent diagonalization of the HF Hamiltonian. Obviously that for the periodic gauge 
SCk. M qo ,<,,, ko iS2|S3 ) = 1, and this formulation is automatically reduced to the standard one 
(Ref. 3). The correction (j7|) was not considered in the previous works (Ref. and 8). 
However, it is found to be small (basically, within the accuracy of the model analysis). Nev- 
ertheless, in a number of cases, the extension ([7]) can be important in order to reproduce 
the correct symmetry of the vector P. 

Moreover, we have found a numerical error in the previous calculations of P, reported 
m Ref. I and 3 all numerical values reported in these two papers should be additionally 
divided by about 2.5. This error was corrected in the present calculations. 

Finally, as will become clear below, the main contribution to the electronic part of AP 
is magnetic and is related to the formation of some complex magnetic texture. Therefore, 
we will also call it 'the magnetic part of AP'. Moreover, as we will see in Sec. 1111} the ionic 
polarization, parallel to the orthorhombic a axis, is associated with the oxygen and yttrium 
sites. Thus, one can clearly separate the contributions of different atomic site: yttrium and 
oxygen atoms contribute to the ionic part of AP a , while Mn atoms - to the magnetic part of 
AP a (thereafter, we will denote by AP a the electric polarization parallel to the orthorhombic 
a axis). 



B. Behavior of interatomic magnetic interactions 



In this section we briefly discuss the effect of the FE displacement on the behavior of 
isotropic exchange interactions, without the relativistic SO coupling. 

Parameters of interatomic magnetic interactions, derived using the perturbation the- 
ory expansion for infinitesimal spin rotations near the FM state,— ^ are listed in Ta- 
ble HI These parameters correspond to the following definition of the spin Hamiltonian: 
Hs = — Yl/ij) Jij e i ' e ji where ej is the direction of the spin at the site i, and the summation 
runs over the inequivalent pairs (ij). In agreement with the previous finding,-^ the magnetic 
structure in the ab plane is mainly controlled by Ji, J^, and J3. 

The FE displacements in the P2\nm structure mainly affect the nearest-neighbor (NN) 
interactions J\. Namely, since FE displacements create two types of the Mn-Mn bonds, 
alternating along the b axis (see Fig. [1]), the corresponding magnetic interactions become 
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TABLE I. The values of interatomic magnetic interactions (in meV) in the ab plane of YMnC>3, 
obtained for the experimental P2\nm and Pbnm structures. The notations of magnetic interactions 
are explained in Fig. [TJ 



crystal structure 


J? 


J? 


T a 

J 2 




h 


P2\ nm 


-5.8 


-3.7 


-0.2 


-1.2 


-2.2 


Pbnm 


-5.0 


-5.0 


-0.2 


-1.2 


-2.3 



also inequivalent and J±, operating in the compressed bond, is generally different from 
J{" , operating in the stretched one. Then, we have obtained that these interactions satisfy 
the following inequality: \J±\ > \Ji\, meaning that the compression of the Mn-Mn bonds 
strengthens the AFM coupling (and vice versa). Such a behavior can be rationalized by 
considering the double exchange (DE) and superexchange (SE) contributions to J\. Of 
course, the concept of the DE physics is rather unusual for insulators. We understand it in 
the following sense: generally, in order to find the energy gain, caused virtual hoppings in an 
insulating state, one should employ the SE theory. Then, in the ferromagnetically coupled 
bond, this energy gain is proportional to t 2 /A (where t is an effective transfer integral, and 
A is the splitting between occupied and empty levels with the same spin). However, in band 
insulators, A is created by the same hybridization effects as t. Therefore, A ~ t, and the 
energy gain is also proportional to t. In this sense, the mechanism can be called the 'double 
exchange'. The semiquantitative theory of DE and SE interactions can be formulated in 
terms of the 1/A ex expansion for Ji,— where A cx is the intraatomic exchange splitting 
between the majority- and minority-spin states (note that in manganites, A ex is the largest 
parameter: ~4.7 eV in the case of YMnOs). Then, the DE describes the zeroth-order effects 
with respect to 1/A CX and the SE - first-order effects. Alternatively, one can say that the 
DE and SE are proportional to transfer integrals, in the first and second order, respectively. 



The details can be found in Ref. 
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Following this strategy, the DE interactions in the P2\nm structure can be evaluated as 
23.9 meV and 23.2 meV, for the compressed and stretched bonds, respectively (23.8 meV in 
the Pbnm phase). The SE interactions are of the order of —30.7 meV and —27.9 meV, for 
the compressed and stretched bonds, respectively (—29.6 meV in the Pbnm phase). Thus, 
the change of Ji, as well as difference between J< and J{" in the FE state is mainly controlled 



by the superexchange processes. 

Similar tendency was found for the NN interactions between adjacent ab planes, where 
the compression of the bonds additionally stabilizes the AFM coupling. For example, the 
interlayer coupling, operating in the compressed and stretched bonds of the P2\nm phase, 
was found to be —7.3 meV and —6.8 meV, respectively (in comparison with —7.1 meV, 
obtained for the Pbnm phase). 

C. Ferroelectric displacements, stability of the E-state, and magnetic contribution 
to the electric polarization 

Fig. [2] shows results of HF calculations for the general twofold periodic magnetic texture 
without SO coupling. In this setting, the spins in each magnetic sublattice are aligned 
antiferromagnetically along the b axis, and two such sublattices in the ab plane are allowed 
to rotate relative to each other (meanwhile, the magnetic coupling along the c axis was 
fixed to be AFM). The relative direction of spins in two magnetic sublattices is specified 
by the angle <fi. Then, <fi = 90° corresponds to the spin-spiral alignment, while = and 
180° correspond to two AFM configurations of the E-type, which are degenerate in the 
centrosymmetric Pbnm structure, as was discussed in Ref. y. The details of calculations can 
be found in Ref. [fl 

As expected, the FE displacements lift the degeneracy of the E-states and stabilize only 
one of them (corresponding to <p — 180° in our setting). In the ground-state configuration, 
the AFM and FM coupling is set in those Mn-Mn bonds, which are compressed and stretched, 
respectively. This is consistent with the behavior of interatomic magnetic interactions, 
discussed in Sec. Ill BL The total energy difference between stable (<f) = 180°) and unstable 
(<ft = 0) E-states is about 4.1 meV/Mn. For the Pbnm structure, the total energy minimum 
corresponds to the spin spiral configuration, while the E-states are higher in energy by about 
0.08 meV/Mn.- Thus, the FE displacements have a profound effect on the total energy of 
orthorhombic YMnC>3. 

On the contrary, the magnetic state dependence of electronic polarization is rather insen- 
sitive to whether the FE displacements are present or not. In the centrosymmetric Pbnm 
phase, the electronic polarization identically vanishes for the spin-spiral state, AP|' = 0, 
as expected without the SO coupling.- In the P2^nm phase, AP.^(90 o ) is finite. However, 
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FIG. 2. (Color online) Results of Hartree-Fock calculations for the low-energy model, constructed 
using experimental parameters for the P2\nm and Pbnm structure. Upper panel shows the be- 
havior of electronic polarization for the general twofold periodic magnetic texture, depending on 
the angle (j) between Mn-spins in two magnetic sublattices. The polarization is measured relative 
to the spin-spiral state. The magnetic structure in the ab plane is explained in the inset. Lower 
panel shows corresponding dependence of the total energy on the angle <j>. The magnified behavior 
of the total energy for the Pbnm structure is shown in the inset (the notations and the scale of the 
axes are the same as in the main figure). 

its absolute value is more than one order of magnitude smaller than in the E-state. Thus, 
to a good approximation, AP.^(90 o ) can be used as the reference point in the analysis of 
the magnetic-state dependence of AP^'(0), both in the Pbnm and P2inm phase. This 
magnetic-state dependence is shown in Fig. [2]^a). One can clearly see that the behavior of 
AP^'(0) is practically identical for the Pl^nm and Pbnm structures. 

Thus, one can expect that the FE displacements in the P2 1 nm phase may have only 
indirect effect on the magnetic part of the electric polarization: they define the equilibrium 
value of <p, which controls AP^ 1 . The direct contribution of the FE displacements to AP^ 1 is 
small. Nevertheless, they will contribute to the ionic polarization, which will be discussed 
below. 
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FIG. 3. (Color online) Distribution of spins in the ab plane of twofold periodic magnetic texture 
of YMnC-3, as obtained for the experimental Pbnm (a) and P2\nm (b) structure. The manganese 
atoms are indicated by the big red (dark) spheres and the oxygen atoms are indicated by the small 
green (grey) spheres. The numbers stand for the angles between spin magnetic moments in the 
Mn-O-Mn bonds. 



D. Magnetic ground state 



In this section, we investigate the effect of the relativistic SO coupling on the magnetic 
ground state of YMnCV We start from the collinear E-state, switch on the SO interaction, 
and further optimize the magnetic texture, by solving the electronic low-energy model ([T]) in 
the HF approximation (the details can be found in Ref. |6|). The results of this optimization, 
which was performed for the experimental P2 1 nm and Pbnm structures, are explained in 
Fig. |3j In both cases, the SO coupling leads to a noncollinear alignment. Nevertheless, 
apart from some quantitative differences, two magnetic textures look rather similar. Thus, 
the magnetic interactions, operating in the Pbnm structure, are already sufficient to form 
the noncollinear magnetic texture, which breaks the inversion symmetry. The relativistic SO 
interaction plays a crucial role in this inversion symmetry breaking. Indeed, let us consider 
the spin texture, which was obtained for the Pbnm structure (Fig. [3^). The spins in each 
magnetic sublattice are ordered antiferromagnetically along the b axis. Such a magnetic 
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order is stabilized by the magnetic interactions and J3 (see Fig. [T]). 9 Then, in each 
magnetic sublattice, the spins tend to align parallel to the longest Mn-0 bonds and thus 
minimize the single-ion anisotropy energy. In the combination with the AFM coupling along 
the b axis, this will naturally lead to the alternation of angles between neighboring spins 
(60° and 120°, respectively) and the formation of the two inequivalent Mn-O-Mn bonds. As 
soon as the bonds become inequivalent, the lattice will relax in order to adjust this change of 
the symmetry. This qualitatively explains the origin of the FE displacements in the P2\nm 
structure. Of course, the change of the crystal structure from Pbnm to P2inm will further 
affect the magnetic texture, which becomes 'more collinear' (the angles between neighboring 
spins are 20° and 160°, respectively). 

Another important issue is the competition between twofold and fourfold periodic mag- 
netic textures. The (nearly) fourfold periodic magnetic structure is realized, for example, in 
TbMnCV- In the previous work, one of us has argued that, in the centrosymmetric Pbnm 
phase, the fourfold periodic magnetic texture tends to have lower energy (both for TbMn03 
and HoMnOs, although experimentally these two materials are known to form fourfold and 
twofold periodic textures, respectively).- Therefore, it is reasonable to ask how this tendency 
will be affected by the FE displacements in the Pl^nm phase. Intuitively, we expected that 
the FE displacements should stabilize the twofold periodic magnetic texture, in agreement 
with the experiment. 

The fourfold periodic magnetic solution was obtained by starting from the spin-spiral 
configuration (the so-called ab helix). Then, after switching on the SO coupling, the spin 
spiral steadily equilibrates to the magnetic textures, which are shown in Fig. HI We would 
like to emphasize that there is a substantial deviation from the spin-spiral alignment both for 
the Pbnm and P1\nm phases. The obtained textures are inhomogeneous in a sense that, 
like in the twofold periodic case, they are also characterized by the alternation of angles 
between neighboring spins in the ab plane. For example, the magnetic textures in the 
P2\nm phase is characterized by three such angles: 128°, 129°, and 156° (for comparison, 
in the homogeneous spin spiral, all angles are equal to 90°). As was argued in Ref. |&, this 
inhomogeneity plays a very important role and is responsible for the FE activity in the 
fourfold periodic magnanites, similar to the twofold periodic ones. Thus, we believe that 
there is no conceptual difference between twofold and fourfold periodic systems: 

• both superstructures have a magnetic origin, where the relativistic SO interaction 
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FIG. 4. (Color online) Distribution of spins in the ab plane of fourfold periodic magnetic texture 
of YMnOs, as obtained for the experimental Pbnm (a) and P2\nm (b) structure. The manganese 
atoms are indicated by the big red (dark) spheres and the oxygen atoms are indicated by the small 
green (grey) spheres. The numbers stand for the angles between spin magnetic moments in the 
Mn-O-Mn bonds. 

plays an important role; 

• in both superstructures, the FE activity is caused by the inhomogeneity of the magnetic 
texture or by its deviation from the conventional spin-spiral alignment. 

In the Pbnm phase of YMn03, the fourfold periodic texture has lower energy than the 
twofold periodic one. This is consistent with results of the previous work.- The total en- 
ergy difference between these two states, obtained in the HF approximation, is about 2.6 
meV/Mn. Somewhat unexpectedly, we have found that, in the P2inm phase, the fourfold 
periodic texture still has lower energy, although the total energy difference is reduced from 
2.6 till 1.0 meV/Mn. Moreover, the distribution of magnetic moments, obtained for these 
two crystallographic phases, is very similar (see Fig. H]). Thus, like in the twofold periodic 
case, there is only a quantitative difference between fourfold periodic magnetic textures, 
realized in the Pbnm and P2\nm phases. This difference is also quite modest: for example, 
the angles between neighboring spins change only within 10-20°. 
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Probably, the precise analysis of the delicate balance between different magnetic states 
is beyond abilities of our low-energy model, because the model itself does not includes 
several ingredients, as was discussed in Ref. |9|. Nevertheless, it clearly shows that the FE 
displacements in the P2inm phase work in the 'right direction' and tend to additionally 
stabilize the twofold periodic magnetic texture, in agreement with the experiment. Moreover, 
the fact that the twofold and fourfold periodic magnetic states are so close in energy may lead 
to the coexistence of these states, as was proposed recently in some experimental works.— 



E. Values of ferroelectric polarization 

In this section we analyze different contributions to the FE polarization in the Pbnm and 
P2±nm phase of YMnC>3. Obviously, the ionic contribution vanishes in the centrosymmetric 
Pbnm phase, but becomes finite in the Pl^nm one. Therefore, in the P2\nm structure, the 
definition of the electronic and ionic terms should be consistent with each other and also 
with the direction of the FE displacement. Then, the first question is which values of the 
atomic charges {Z{\ shall we use in our model analysis? The first possibility is to take the 
ionic charges: i.e., — 2\e\ for O and +3|e| for Y and Mn. In this case, the ionic contribution 
to AP a can be estimated as 0.55 /zC/cm 2 . The second possibility is to take {Z{\ from the 
band structure calculations in the linear muffin-tin-orbital (LMTO) method,— where each 
Zi is understood as the the positive nuclear charge plus negative electronic charge within the 
atomic sphere i (including empty spheres, which are used for a better space filling). Then, 
the first term in Eq. fl3]) can be estimated as —0.04 /^C/cm 2 . This simple analysis shows that 
the ionic part of AP a is sensitive to the definition of {Zi}. Anyway, it seems to be smaller 
than the electronic contribution to AP a , caused by the noncentrosymmetric distribution of 
spins. 

The electronic polarization, obtained for the twofold periodic magnetic texture in the 
P2\nm phase, is about —0.91 /zC/cm 2 . Then, the total polarization AP a can be estimated 
as —0.95 -I — 0.36 /iC/cm 2 , depending on the approximation, employed for the ionic part 
of AP a . These values are somewhat overestimated in comparison with the experimental 



data: | AP a | = 0.24 /xC/cm 2 in the bulk (Ref. |7|) and 0.80 /xC/cm 2 in the thin films (Ref. Il6j). 
Nevertheless, the disagreement is quite modest: at least, we do not encounter the notorious 
'one order of magnitude difference', which takes place if one uses the theoretically optimized 

15 



crystal structure instead of the experimental one.->£ We will come back to the analysis of 
this problem in Sec. IIHI Moreover, one should keep in mind two points: 

(i) The electronic polarization in the twofold periodic texture is very sensitive to the 
angle <fi between spins in two magnetic sublattices. Moreover, the closer the system to 
the spin-spiral alignment is, the smaller the value of |AP a '| will be (Fig. [2]). For example, 
in the Pbnm phase, <fi is closer to 90° (Fig. [3]). Therefore, the electronic polarization is 
substantially smaller: |AP.^|= 0.55 /iC/cm 2 . Nevertheless, due to the FE displacements 
in the P2\nm phase, the magnetic alignment becomes more collinear and the value of the 
electronic polarization increases. From this point of view, one reason why the experimental 
value of |AP a | is somewhat smaller than the theoretical one may be related to the fact 
that the experimental magnetic texture itself is more noncollinear. It would be interesting 
to check it experimentally, by measuring the deviation from the collinear E-type AFM 
alignment. 

(ii) Another mechanism, which might reconcile disagreement with the experimental data 
is the coexistence of the two- and fourfold periodic magnetic textures, which was suggested 
recently for the thin films of YMnC>3.— According to our model analysis, the realization of 
the fourfold periodic magnetic texture in the P1\nm phase (instead of the twofold periodic 
one) would lead to the drop of the electronic polarization from —0.91 till 0.01 /iC/cm 2 . 
Thus, the mixture of these two magnetic states will probably decrease |AP a |. 



III. STRUCTURAL OPTIMIZATION: RESULTS AND ANALYSIS 

In this section we report results of structural optimization for the orthorhombic YMnC>3, 
by assuming the collinear E-type AFM alignment, and discuss possible implications of the 
obtained structure on the magnetic and FE properties. Namely, our strategy is the following. 
First, we perform the structural optimization by using both LSDA and LSDA+£/ functionals 
without the relativistic SO coupling. The main purpose of these calculations is to investigate 
the effect of the exchange striction, caused by the collinear E-type AFM alignment. In this 
study, the value of U was treated as a parameter. Namely, we considered two such values: 
U— 2.2 eV - an averaged value, which is typically obtained for the low-energy model of 
orthorhombic manganites,- 1 ^ and U= 6.0 eV, imitating the 'large-?/ limit' (at least, the FE 
polarization in this regime is nearly saturated - see Refs. OandH). As for the intraatomic 

16 



exchange coupling, we take the atomic value Jh= 0.9 eV. 

For the structural optimization, we use the projector augmented-wave method imple- 
mented in the VASP code package.— The k-space integration was performed on the mesh of 
the 6x6x6 points in the first Brillouin zone. The plane wave cut-off energy was chosen to be 
500 eV. The optimization was performed by starting from the centrosymmetric experimental 
Pbnm structure. 

Then, we take the optimized crystal structure and again construct the low-energy model 
for each set of the parameters, obtained for different values of U. After that, we solve the 
low-energy model and discuss the relative stability of the two- and fourfold periodic textures 
as well as the different contributions to the FE polarization. Thus, we use LSDA+?7 basically 
as a 'machinery', which produces some crystal structure, and then we want to answer the 
question how good is this crystal structure, by analyzing more delicate magnetic properties 
and comparing them with the ones obtained for the experimental structure. 

In Table [TTJ, we show results of structural optimization with U— 2.2 eV. Similar results, 
obtained in LSDA and LSDA+£7 with U= 6.0 eV, are collected in Ref. [ll|. At first sight, 

TABLE II. Parameters of the P2\nm structure, obtained in LSDA+C/ with U= 2.2 eV for the E- 
type antiferromagnetic ordering. In this table, x, y, and z are the fractional coordinates, calculated 
in terms of the lattice parameters a= 5.182 A, b= 5.740 A, and c= 7.251 A, respectively. 



site x y z 

Yl 0.478 0.831 1/2 

Y2 0.520 0.663 

Mn 0.755 0.251 

01 0.294 0.569 0.303 

02 0.887 0.717 1/2 

03 0.110 0.791 

04 0.698 0.927 0.196 



the internal coordinates, measured in units of the orthorhombic lattice parameters, seem to 
agree with the experimental values reasonably well.- 1 ^ A more serious discrepancy occurs 
for the unit-cell volume, which is underestimated by about 6 % in LSDA and by 2 % in 
LSDA+[7 with U— 6.0 eV. The values of Mn-Mn bondlengths and Mn-O-Mn angles in the 

17 



ab plane, obtained for different values of U, are summarized in Table IHIl The notations 
are explained in Fig. |5j Although there is a correlation between the theory and the exper- 
iment (for example, the Mn-Mn bondlengths and Mn-O-Mn angles are generally larger for 
the FM bond, etc.), none of the calculations reproduces the experimental values. LSDA 
overestimates the FE displacements: both Mn-O-Mn angle and Mn-Mn bondlength are sub- 
stantially smaller for the the AFM bond in comparison with the FM one. This problem 
can be resolved by including the Coulomb U, but only partially. For example, the best esti- 
mate for the AFM Mn-Mn bondlength was found for U= 6.0 eV. However, the FM Mn-Mn 
bondlength practically does not depend on U. Both values are smaller than the experimental 
ones, presumably because the unit-cell volume is also smaller. Moreover, the calculations 
overestimate the difference between the Mn-O-Mn angles in the FM and AFM bonds, even 
for U= 6.0 eV. 

TABLE III. Results of structural optimization for the E-type AFM structure in LSDA and 
LSDA+[/ with U = 2.2 and 6.0 eV: the Mn-Mn distances (d, in A) and Mn-O-Mn angles (Z, 
in degrees), as obtained for the FM and AFM bonds in the ab plane of YMnOs, in comparison 
with the experimental values (exp) reported in Ref. [7J. The geometry of the bonds is schematically 
explained in Fig. [5j 







AFM bond 




FM bond 




d(Mn-Mn) 


Z (Mn-O-Mn) 


d(Mn-Mn) 


Z (Mn-O-Mn) 


LSDA 


3.75 


143.8 


3.88 


148.5 


U= 2.2 eV 


3.82 


143.7 


3.91 


146.8 


U= 6.0 eV 


3.87 


143.9 


3.89 


145.2 


exp 


3.91 


143.8 


3.94 


144.1 



Thus, it is interesting to ask how all these changes of the crystal structure are reflected 
in the behavior of the FE polarization and the magnetic properties of YMn03? 

The structural optimization in LSDA and LSDA+£7 provides a reasonable starting point 
for the analysis of the noncollinear magnetic texture. At least, the results are consistent with 
each other and can be interpreted in the following way. Since LSDA overestimates the FE 
displacements, the obtained twofold periodic magnetic texture is nearly collinear: the angles 
between NN spins in the ab plane are 6° and 174° (Table HV|) . which are close to 0° and 180°, 
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FIG. 5. (Color online) Schematic representation of the Mn-O-Mn bonds in the E-type AFM 
structure. The numerical values for the Mn-Mn distances and Mn-O-Mn angles denote results of 
structure optimization in LSDA+?7 with U= 2.2 eV. 

realized in the collinear E-state. The Coulomb U decreases the FE displacements. Therefore, 
the noncollinearity will systematically increase, and U= 6.0 eV practically reproduces the 
values of 0, obtained for the experimental P2\nm structure. A similar behavior was found 
for the fourfold periodic magnetic texture: the LSDA parameters of the crystal structure 
lead to a nearly collinear magnetic alignment (the angles between NN spins are 11°, 176°, 
and 178°, respectively - Table HV|) . while the noncollinearity again systematically increases 
with the increase of U. 

However, the situation with the FE polarization appears to be more serious. Let us start 
with the analysis of the ionic contribution, which was computed by using ionic charges. 
Table El shows partial contributions to this polarization, associated with different atomic 
sites. Then, we can readily answer the question why the polarization, obtained for the 
experimental P2\nm structure, is relatively small (only 0.55 /iC/cm 2 ). This is because 
there is a strong cancelation of contribution of the oxygen and yttrium sites (note, that the 
contributions of manganese sites to the ionic polarization, parallel to the orthorhombic a 

19 



TABLE IV. Results of solution of four low-energy models, constructed for the experimental P2\nm 
structure (exp) and theoretically optimized structures in LSDA and LSDA+J7 with U= 2.2 eV 
and U= 6.0 eV. The structural optimization was performed by imposing the collinear E-type 
antiferromagnetic alignment. In this table, <j> denotes the angles between neighboring Mn-spins in 
the ab plane of twofold and fourfold periodic magnetic texture (denoted as '2-fold' and '4-fold', 
respectively), AP^ 1 is the value of electronic polarization parallel to the orthorhombic a axis, and 
AE is the total energy difference between two- and fourfold periodic magnetic textures. For the 
experimental P2inm structure, the distribution of the angles eft is explained in Figs. 0and[H 



property 


exp 


LSDA U= 2.2 eV 


U= 6.0 eV 


2-fold <p (°) 


20, 160 


6, 174 8, 172 


23, 157 


2-fold AP? (/jC/cm 2 


) -0.91 


1.66 0.86 


1.18 


4-fold <p (°) 


128, 129, 156 


11, 176, 178 31, 167, 171 


68, 157, 157 


4-fold AP? (fiC/cm 2 


) 0.01 


0.77 0.36 


0.43 


AE (meV/Mn) 


1.0 


-3.2 -0.4 





TABLE V. Total and partial contributions 


to the ionic polarization parallel to the orthorhombic a 


axis, originating from the displacements of 4 yttrium atoms (4Y), 8 oxygen atoms in the ab plane 


(80 a b, denoted as 01 and 04 in Table ILT|) . and 4 oxygen atoms between the planes 


(40 c , denoted 


as 02 and 03 in Table HI|) . which were obtained for the experimental P2\nm structure (exp) and 


theoretically optimized structures in LSDA and LSDA+C/ with U= 2.2 eV and U= 6.0 eV. All 


values are in /xC/cm 2 . 


The structural optimization was performed by imposing the collinear E-type 


antiferromagnetic alig 


nment. 






contribution 


exp 


LSDA U= 2.2 eV 


U= 6.0 eV 


4Y 


-0.35 


0.14 -0.53 


0.08 


80 ab 


0.73 


4.34 2.61 


1.15 


40 c 


0.17 


0.09 0.33 


0.03 


total 


0.55 


4.57 2.41 


1.26 



axis, vanish, because all Arj||a =0). In the theoretical calculations, a similar situation was 
found only for U— 2.2 eV, while in LSDA and for U— 6.0 eV all contributions are positive. 
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Moreover, the contributions of the oxygen atoms in the ab plane are strongly overestimated. 
Therefore, the ionic polarization, obtained for the optimized theoretical structures, is already 
larger than the experiment one. 

The disagreement becomes even more serious when take into account the electronic polar- 
ization (Table HV|) . which for all theoretical structures is aligned in the same direction as the 
ionic one. A similar situation was encountered in the previous calculations for HoMnC>3.- 
This behavior is clearly different from the one obtained for the experimental P2\nm struc- 
ture, where the ionic and electronic terms were antiparallel. 

The additional noncollinearity, obtained for the less distorted structures (in the regime of 
large U), decreases AP|'. This tendency was found both for the two- and fourfold periodic 
magnetic texture (Table II Vj) . Nevertheless, this behavior is nonmonotonous, and there is 
something special in the theoretical L U= 2.2 eV structure. First, only for this theoretical 
structure, there is a partial cancelation of ionic contributions, associated with the oxygen 
and yttrium atoms, similar to the experimental structure. Then, this theoretical structure 
produces the minimal value of AP® 1 , which is also favorable from the viewpoint of com- 
parison with the experimental polarization. One the basis of these observations, one may 
speculate that U= 2.2 eV generates the best theoretical crystal structure. However, the 
total polarization is largely overestimated, and the agreement with the experimental data is 
still poor. 

The use of the LMTO charges, instead of the ionic ones, does not resolve the problem: 
although the ionic polarization becomes somewhat smaller (about 0.35 /iC/cm 2 for U= 6.0 
eV), it has the same sign as AP^ 1 . Thus the cancelation does not occur. Moreover, in the 
ionic polarization, the contributions of the yttrium and oxygen atoms also have the same 
sign and do not cancel each other. 

Apparently, the partial cancelation of the electronic and ionic polarizations, as well as 
the contributions of the oxygen and yttrium sites to the ionic part, explains relatively small 
value of AP a , which was observed in the experiment. It also indicates at a serious problem, 
existing in the electronic structure calculations, which fail to reproduce some fine details 
of the crystal structure, responsible for this cancelation. One reason may be related with 
fundamental problems of LSDA or LSDA+[7 (for instance, the problem of incorrect charge- 
transfer energy in LSDA+[7, which cannot be resolved simply by changing the value of U 
- Ref. Q). Another possibility is that the structural optimization should be performed 
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for the noncollinear magnetic texture, after including the SO coupling. However, from the 
computational point of view, this procedure is extremely difficult. For example, even without 
the structural relaxation, the optimization of the magnetic structure in the low-energy model 
requires several thousands of iterations. Obviously that simultaneous optimization of both 
crystal and magnetic structure in the full scale electronic structure calculations will be much 
more computationally demanding. 



IV. SUMMARY AND CONCLUSIONS 



We discussed the origin of the multiferroic state in the twofold periodic orthorhombic 
manganites, by considering YM11O3 as an example. Traditionally, the FE activity in these 
systems is attributed to the exchange striction effects, associated with the collinear E-type 
AFM alignment. In this picture, the FE displacements lift the frustration of isotropic ex- 
change interactions and stabilize the E-type AFM state, while the relativistic SO interaction 
is believed to be irrelevant to this type of processes. 

Instead, we argued that the E-state can be regarded as a combined effect of frustrated 
isotropic interactions and the single-ion anisotropy energy, acting in the ab plane of mangan- 
ites. The SO interaction plays a crucial role in this case: it lifts the degeneracy and stabilizes 
a canted E-type AFM state. Formally speaking, in order to form the noncentrosymmetric 
E-type AFM state, it is sufficient to have the centrosymmetric Pbnm structure, while the FE 
displacements in the P1\nm structure only follow this lowering of the magnetic symmetry. 
As far as the equilibrium magnetic structure is concerned, the use of the P2\iim structure 
instead of the Pbnm one yields only modest quantitative changes, while all main trends can 
be obtained already in the centrosymmetric Pbnm phase. 

In this analysis, we employed the realistic low-energy model,— derived from the first- 
principles electronic structure calculations,— and the experimental parameters of the crystal 
structure.- 

Then, we addressed an important issue of the 'order of magnitude difference' between the 
theoretical and experimental values of the FE polarization AP a , which is typically observed 
for the twofold periodic systems. 3 We have found that, if one uses the experimental P2\nm 
structure, one can get a reasonable agreement also with the experimental value for the FE 
polarization (similar conclusion was drawn recently in Ref. |7|). The reason for it is the partial 
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cancelation of different contributions in AP a , which takes place for the experimental P2\nm 
structure. Particularly, there is a cancelation of electronic and ionic parts, as well as the 
contributions of the oxygen and yttrium atoms to the ionic polarization. 

In the last part of this work, we attempted to explore the same questions on a purely 
theoretical level, and derived parameters of the noncentrosymmetric P2inm phase from 
the first-principles calculations. For these purposes we employed the LSDA and LSDA+[7 
techniques without the relativistic SO coupling, and performed the structural optimization 
by assuming the collinear E-type AFM alignment. After that, we used again the obtained 
structural information for the construction and solution of the low-energy model. We have 
found that, although the obtained theoretical structure stabilizes the canted E-type AFM 
state, it fails to reproduce the experimental values of the FE polarization. Thus, it seems 
that whenever we use the theoretical parameters of the crystal structure, we face the 'order 
of magnitude difference' problem for the FE polarization. Namely, the cancelation, which 
took place for the experimental P2\nm structure, does not occur for the theoretical one. 
Intuitively, this result is in line with our main idea: if the SO interactions plays a crucial role 
in stabilizing the (canted) twofold periodic E-state, it cannot be ignored in the process of 
the structural optimization. On the other hand, there may be more fundamental problems, 
related to the LSDA and LSDA+[7 functionals. We hope that our analysis will stimulate 
further theoretical and experimental works aiming to resolve this interesting and important 
issue. 
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